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Abstract 

The main objective of this paper is to show that two asymptotically 
stable steady states which belong to an analytic path of asymptotically stable 
steady states can be gradually transferred one to the other by successive 
changes of the control parameters. 

1 Introduction 

For nonlinear systems of differential equations with control, it has been proved 
(see P), that two asymptotically stable steady states belonging to an analytic 
path of asymptotically stable steady states can be transferred one in the other 
by successive maneuvers along the path. That means, according to [llj . that the 
process described by the system can be piloted through the domains of attraction 
of the intermediary steady states, from the first to the second steady state. 

In this paper, a similar result is established for the nonlinear discrete dynamical 
systems with control. A theorem from 4^ is used, which states that the domain of 
attraction of an asymptotically stable fixed point of a nonlinear system of analytic 
difference equations, is the natural domain of analyticity of a certain Lyapunov 
function. 



2 Preliminaries 

We consider the following nonlinear discrete dynamical system with control: 

Xk+i^f{xk,a) fc = 0,l,2... (1) 

In CJ, / is a given function f : Q x D ^ , Q d , D d M™ are domains, 
a; € n is the state parameter, and a E D is the control parameter. What concerns 
the regularity of /, we assume that / is an analytic function. 
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A state a;° G is a steady state for if there exists a^o s D such that 

x° = /(x°,a,o) (2) 

The steady state G f2 of is stable''' provided that given any ball 
B(x^,e) = {x€ n/\\x-x°\\ < e}, there is a ball B{x°,S) ^ {x £ n/\\x-~x^\\ < S} 
such that if a; e B{x°,S) then f'^{x,axo) € B{x°,e), for fc = 0, 1,2, ... ^5^. 

If in addition there is a ball r) such that f'^{x, a.j;o) ^ as fc oo for 

all X € B{x^,r) then the steady state x^ is ''^asymptotically stable'^ [S]. 

The domain of attraction DA{x^) of the asymptotically stable steady state x'^ 
is the set of initial states x G from which the system converges to the steady 
state itself i.e. 

DA[x°) = {xe n\!''{x, a,o) x"} (3) 

An analytic path of steady states of (^) is an analytic function : Di C I? — » O 
which satisfies 



tp{a) = f{ip{a),a), for any a £ Di (4) 

An analytic path of asymptotically stable steady states of is an analytic path 
of steady states which are all asymptotically stable. 

A change of control parameters from a' to a" in |^ is called maneuver and is 
denoted a' — > a". The maneuver a' a" is successful on the path ip if a' , a" G Di 
and the sequence defined by 

Xk+i = fixk,a"), xo = ip{a') (5) 

tends to ip{a") as cx). 

The following proposition from 0] concerning the discrete dynamical systems 
without control parameters is helpful. 

Proposition 1. // the analytic function 5 : A — > A from the system 

yk+i=9(yk), /c = l,2,... (6) 

satisfies the following conditions: 

5(0) = (7) 



\\do9\\ < 1 (8) 

then is an asymptotically stable steady state of DA(0) is an open subset of 
A and coincides with the natural domain of analyticity of the unique solution V 
of the iterative first order functional equation 

V{g{y))-V{y)^-\\yr 
F(0) = 



The function V is positive on DA{0) and V{y) — > +00, for any G FrDA{0) 
(FrDA{0) denotes the boundary of DA{0)). 
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Remark 1. If tp : Di C D ^ Q. is an analytical path of steady states of Qj, then 
the function g(-, a) : Q — >f{a) fl — (p{a), given by 

9{y,a) ^ f{y + (p{a),a) ~ ip{a) for y e n - ip{a) (10) 

is analytic and satisfies g{0, a) — 0, for any a G Di. Therefore, y = is a steady 
state of the system 

yk+i=giyk,a) fc = 0, 1,2... foranyaGDi. (11) 

The steady state x ~ <p(a) of 0) is asymptotically stable if and only if the 
steady state y = of the system is asymptotically stable. The relationship 
between the domain of attraction of the steady state x — (p{a) of Q) and that of 
the steady state y ^ of |77J) is DA{ip{a)) = ip{a) + DA{0). 

3 Theoretical results 

We now state an existence theorem for an analytic path of asymptoticaUy stable 
steady states of (QJ. 

Theorem 1. If the analytic function f from 0) satisfies: 

1. there exist (a;°, a^) e fl x D such that — /(a;°, a°) 

2. ||a,o/(xO,aO)|| < 1 

then there exists a maximal domain Di C D containing and a unique ana- 
lytic path ip : Di —^Q,of asymptotically stable steady states of 0j satisfying the 
following conditions: 

a. (^(a") = x*^ ; 

b. \\dxf{p{a),a)\\ < 1 for any a G Di; 

c. For a' , a" G Di the maneuver a' — > a" is successful on the branch ip if and 
only ifip{a') belongs to the domain of attraction of Lp{a"). 

Proof. As the functions / and dxf are continuous on x D, taking into account 
the properties 1. and 2. of /, there exist two maximal domains i7i C and Di 
and a unique analytic function p : Di — + fli such that: 

1. {x°,a°) eflix Di and ip{a") = 2;°; 

2. Lp{a) = f{(p{a),a), for any a G -Di; 

3. \\da:f{(p{a),a)\\ < 1 for any a G Di. 

This means that ip is path of asymptotically stable steady states for ^ (see 
Proposition 1 and Remark 1). 

A maneuver a' —^ a" is successful on the path p if and only if the sequence 
given by © tends to p{a") as fc ^ 00, which means that p{a') belongs to the 
domain of attraction of p{a"). □ 
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In the followings, it is assumed that the conditions of Theorem 1 are fulfilled 
and thus, there exists an analytic path Lp of asymptotically stable steady states of 

Theorem 2. Let he Lp : Di —i-flan analytic path of asymptotically stable steady 
states of 0). There exist an open set G C x D and a non-negative analytic 
function V defined on G satisfying the following conditions: 

a. G D r = {{Lp{a),a)/a £ Di} 

b. 

V{f{x,a),a) - V{x,a) = -||a; - 

V{ip{a),a)=0 ^ ^ 

c. For any a G Di, DA{ip{a)) is the natural domain of analyticity of x —>■ 
V{x, a) 

d. V{x,a) ^-1^5. +00, for any x^ e FrDA{(p{a)) . 

Proof Let be G = IJ iDA{ip{a)) x {a}) C 1^ x L>i and F : G ^ defined by 

qG-Di 



oo 



y(x,a)=^||/'^X2^,«)-^(")f (13) 

Proposition 1 and Remark 1 provide that the set G and the function V{x, a) 
satisfy the conditions a-d. □ 

Corollary 1. If (p : Di fl is an analytic path of asymptotically stable steady 
states of 0) then for any a £ Di there is an open neighborhood Ua of a and an 
open neighborhood U^p(a) of Lp{a) such that: 

1. (p{a') e U^(^a), for any a' G 

2- U,p(a) C DA{Lp{a')), for any a' G Ua 

Proof. For a G Di and x G DA{ip{a)), the function V{x,a) from Theorem 2 
is considered. The real and non-negative function V is defined on the open set 
G = U {DA{(f{a)) X {a}) C 57 x Di, it is continuous and equal to zero on the 

aeDi 

set r = {{ip{a),a)/a G Di} C G. 

As V is continuous and it is equal to zero in {^p(a),a) G G, there is an 
open neighborhood W of (</j(a),a) such that for any (x',a') G W, the inequality 
V{x' , a') < 1 holds. Let be Ua an open neighborhood of a and U^i^a) of fia) such 
that U^(^a) >^ Ua C W. As the function (p is continuous, it can be admitted that 
for any a' G Ua, we have (p{a') G U^(^a) (contrarily, the neighborhood Ua can be 
replaced with a smaller neighborhood Ua C C/q, for which we have ipia') G U^p(^a), 
for any a' G U'a)- 

Thus, for any {x' ,a') G U^(^a) x Ua, we have V{x' ,a') < 1. This means that 
for any x' G ?7ip(Q) and any a' G Ua, we have that x' G DA{p{a')). Thus, 
U^(a) C DA{ip{a')), for any a' e Ua- □ 
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Remark 2. Corollary 1 states that for any a' € Ua, both maneuvers a —f a' and 
a' ^ a are successful on the path ip. 

Theorem 3. For two steady states (p{a*) and ip(a**) belonging to the analytic 
path ip of asymptotically stable steady states of Q), there exist a finite number of 
values of the control parameters , ...,cxP G Di such that all the maneuvers 

a* -^a^ ... -> ^ a** (14) 

are successful on the path (p. 

Proof. Let be P C a polygonal line which joins a* and a**. For any a G P we 
consider the neighborhoods Ua and U^(^a) given by Corollary 1. 

The family of neighborhoods {Ua}a^p is a covering with open sets of the 
compact polygonal line P. From this covering we can subtract a finite covering 

of P, i.e., there exist a^, a^, a* e P such that P C IJ C/gfc. More, it can be 

k=l 

assumed that a* £ C/gi and a** G Uai and that the intersections [/gt flP are open 
and connected sets in P, and 

{Uak n P) n (Uak+2 n P) = for any fc = 1, 2, 17-2. 

Taking into account Remark 2, as a* G J/gi and a** G Uai^ it comes naturally 
that the maneuvers a* and a'^ — > a** are successful on the path ip. 

We still have to prove that each maneuver a'^ — > a*^"*"^ is successful for any 
fc = l,2,...,9-l. 

If a'' G Uak+i, Remark 2 provides that the maneuver a*^ a'^^^ is successful 
on the path ip. 

If a*" ^ Uak+i, a point a'^^'-'+i e [Uak nP)n([/afc+i HP) is considered. Remark 
2 provides that both maneuvers a'' (j^k.k+i ^^^^ ^fe,fc+i _^ q.'=+i are successful 
on the path ip. 

Thus, eventually inserting control parameters a''-'^+^ between a'^ and a'^'^^, 
we come to find (after changing the notation and re-numbering) a finite sequence 
a^jO^, ...,aP € Di such that all the maneuvers 

p -k-k 

a — > a — > a — > ... — > a — > a 

are successful on the path ip. □ 

Remark 3. Theorem 3 states that two steady states belonging to an analytic path 
of asymptotically stable steady states can be transferred one in the other using a 
finite number of successful maneuvers along the considered path. 

4 Numerical examples 
Example 1. 

The following one-dimensional discrete dynamical system with control is consid- 
ered: 

Xk+i = axk{l - Xk), fc = 0,1,2,... (15) 
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where a; £ K is the state parameter and a € M is the control parameter. This 
dynamical system is frequently used for showing chaotic behavior and it is subject 
to several themes of research We will illustrate using this system the concepts 
of analytic path, analytic path of asymptotically stable steady states, domain of 
attraction, successful maneuver. 

For a ^ 0, the steady states of (|15|l arc a; = and x — while for a — 
corresponds only the x — steady state. Thus, for H15() there are three paths of 
steady states: ipi{a) — for a < 0, (p2{a) — 0, for a G K and (/33(a) = 
for a > 0,. The path ipi contains only unstable steady states. The steady states 
belonging to the path (p2 are asymptotically stable for a G (—1, 1). The steady 
states belonging to the path (^3 are asymptotically stable for a G (1, 3). 

In Fig. 11.11 the three paths of steady states are plotted. In Fig. 11.21 the gray 
rectangle represents the reunion of the domains of attraction of the asymptotically 
stable steady states of ips, while the vertical gray line denotes the domain of 
attraction of the steady state (/32(0). In both figures, the black parts of the paths 
of steady states represent the asymptotically stable steady states while the gray 
parts of the paths represent the unstable steady states. 
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Figure 1.1: The paths of steady states for 
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Figure 1.2: The domains of attractions of 
the asymptotically stable steady states of 
the path 1^93 and for the asymptotically 
stable steady state V'2(0) = of (/32 



The domain of attraction of the steady state (p2{0) = is DA{(p2{0)) = while 
the domain of attraction of a steady state <p3(a) for a G (1,3) is DA{(p3{a)) = 
(0,1). These domains of attraction can be obtained using the staircase method 

or can be estimated numerically using the method described in ^). 

The steady state (p3{l.l) can be directly transferred by a single maneuver 
to V53(2.9), because (pail-l) G DA{ip3{2.9)) = (0,1). The DA{ip3{lA)) includes 
1/33(2.9), thus, the maneuver a : 2.9 — > 1.1 is also successful. 

All asymptotically stable steady states of ip^ are in the domain of attraction of 
the steady state f2{0). This means that every maneuver a — + 0, for a G (1,3) is 
successful between the paths tp^ and ip2- Though, a steady state ip2{ct) cannot be 
transferred in an asymptotically stable steady state of (^3, because any maneuver 
of the type a ^ a' , with a' G (1,3) causes a transfer to the unstable steady state 
Ma') = 0. 
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Example 2. 

The following one-dimensional discrete dynamical system with control is consid- 
ered: 

Xk+i = {xk - af + a, A: = 0,1, 2,... (16) 

where a; S M is the state parameter and a e R is the control parameter. 
The sequence Xk, with the starting point xo which satisfies H16|l is: 

Xk ^ {xq - af" + a, /c = 0,l,2,... (17) 

There are three analytic paths of steady states for ((TC|l : (pi{a) — a, ip2{a) = 
a — 1 and ^'3(0^) = a -I- 1, defined for a 6 R. The path ipi is an analytic path of 
asymptotically stable steady states while (^2 and (^3 are analytic paths of unstable 
steady states. In Fig. |31 the continuous line represents the path Lpi, while the 
dashed lines represent the paths 1^92 and (^3. 

For any a G R, the domain of attraction of the asymptotically stable steady 
state •^i{a) is DA{{pi{a)) = [a — l,a+ 1). 

For a* = and a** — 2, let's consider the asymptotically stable steady states 
Lpi{a*) — and ifi{a**) — 2. The maneuver a : a* = — > 2 = a** is not 
successful, because 1^91(0!*) = ^ DA{ipi{a**) = 2) = (1,3). Though, a finite 
number of maneuvers can be found, which transfer the steady state ipi[a*) = to 
the steady state ip\{a**) = 2, for example: 

a : a* = ^ 0.7 ^ 1.4 ^ 2 = a** (18) 
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Figure 2: The paths of steady states for 1)16(1 and the domains of attraction of the 
steady states of </3i corresponding to the control parameters used in ll8l 

These maneuvers are successful, because yi(a*) 6 DA{{pi{Q.T)) = (—0.3, 1.7), 
(pi(0.7) e DA{ipi{lA)) = (0.4,2.4) and (^i(1.4) e DA{ipi{a**)) = (1,2). In Fig. 
the vertical segments represent the domains of attraction of the steady states 
corresponding to the maneuvers ((18|l . 




7 



Example 3. 



The following two-dimensional discrete dynamical system with control is consid- 
ered: 

f Xk+i = {xk ~ a)[{xk - a)^ + ivk - a)"^] + ct ^^g. 
\ Vk+i = iVk - a)[ixk - a)"^ + iVk - a)'^] + a 

where {x, y) G is the state parameter and a £ M is the control parameter. 

There are an infinity of analytic paths of steady states for H19|) : <p(a) = (a, a) 
and ftio) — {a + cost, a + sint), for t G [0, 2tt); all paths are defined for a G M. 
The path ip is an analytic path of asymptotically stable steady states while (pt are 
analytic paths of unstable steady states, for any t G [0, 2tt). 

For any a G K, the domain of attraction of the asymptotically stable steady 
state ip{a) = (a, a) is the ball B{{a, a), 1) = {{x, y) G M^/ [x — aY + [y — aY < !}■ 



Figure 3: The paths of steady states for H19II 

For a* = — 1 and a** = 1, let's consider the asymptotically stable steady states 
f{a*) = (—1,-1) and (p{a**) = (1,1)- The maneuver a : a* = — 1 1 = a** 
is not successful, because ip{a*) = (-1, -1) ^ DA{ip{a**) ^ (1, 1)) = B{{1, 1), 1). 
Though, a finite number of maneuvers can be found, which transfer the steady 
state (p{a*) = (—1, —1) to the steady state Lp{a**) — (1, 1), for example: 

a : a* = -1 ^ -0.5 ^ 0.5 ^ 1 = a** (20) 

These maneuvers are successful, because (p{a*) G DA{ip{— 0.5)) = 0.5, —0.5), 1), 
</j(-0.5) G DA{ip{0)) = B((0,0),1), ip{0) G DA{ip{0.5)) = B((0.5, 0.5), 1), and 
^{0.5) G DAiipia**)) = B{{1,1),1). 



8 



References 



[I] St. Balint: Considerations concerning the maneuvering of some physical sys- 
tems; An. Univ. Timisoara, Ser. St. Mat., vol. XXIII, (1985) p. 8-16. 

[2] M. Fcigcnbaum: Quantitative universality for a class of nonlinear transforma- 
tions; J. Statist. Phys 19, (1978) p. 25-52. 

[3] L. Hormander; An Introduction to Complex Analysis in Several Variables, D. 
Van Nostrand Company, Inc., Princeton, New Jersey 

[4] E. Kaslik, A.M. Balint, S. Birauas, St. Balint: Approximation of the domain 

of attraction of an asymptotically stable fixed point of a first order analyti- 
cal system of difference equations. Nonlinear Studies, I&S Sience Publishers 
(2003). 

[5] W.G. Kelley, A.C. Peterson: Differenc equations; Academic Press, 2001. 

[6] H. Kocak: Differential and Difference Equations Through Computer Experi- 
ments, Springer- Verlag, New York, 1990. 

[7] G. Ladas, C. Qian, P. Vlahos, J. Yan: Stability of solutions of linear nonau- 
tonomous difference equations. Appl. Anal. 41 (1991), no. 1-4, p. 183-191. 

[8] V. Lakshmikantham, D. Trigiante: Theory of Difference Equations: Numerical 
Methods and Aplications, Academic Press, New York, 1988. 

[9] J. LaSalle: Stablity theory for difference equations, MAA Studies in Mathe- 
matics 14, Studies in Ordinary Differntial Equations; edited by J. Hale, 1997, 
p. 1-31. 

[10] J. LaSalle: The Stability and Control of Discrete Processes, Springer- Verlag, 
New York, 1986. 

[II] V. Rasvan: Teoria Stabilitatii, Ed. Stiintifica si Enciclopedica (1987), Bu- 
curesti. 



9 



